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Abstract 

The classical branch-and-bound algorithm for the integer feasibility 
problem 



Finds e QC\Z n , with Q = { x\ 1 < \x < 1 (1) 




has exponential worst case complexity. 

We prove that it is surprisingly efficient on reformulations of (JXJ) , 
in which the columns of the constraint matrix are short, and near 
orthogonal, i.e. a reduced basis of the generated lattice; when the 
entries of A are from {1, . . . , M} for a large enough M, branch-and- 
bound solves almost all reformulated instances at the rootnode. For all 
A matrices we prove an upper bound on the width of the reformulations 
along the last unit vector. 

The analysis builds on the ideas of Furst and Kannan to bound 
the number of integral matrices for which the shortest nonzero vec- 
tors of certain lattices are long, and also uses a bound on the size of 
the branch-and-bound tree based on the norms of the Gram-Schmidt 
vectors of the constraint matrix. 

We explore practical aspects of these results. First, we compute 
numerical values of M which guarantee that 90, and 99 percent of 
the reformulated problems solve at the root: these turn out to be 
surprisingly small when the problem size is moderate. Second, we 
confirm with a computational study that random integer programs 
become easier, as the coefficients grow. 



1 Introduction and main results 



The Integer Programming (IP) feasibility problem asks whether a polyhe- 
dron Q contains an integral point. Branch-and-bound, which we abbreviate 
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as B&B is a classical solution method, first proposed by Land and Doig in 
|20j . It starts with Q as the sole subproblem (node). In a general step, one 
chooses a subproblem Q' , a variable Xi, and creates nodes Q' n {x\xi = 7}, 
where 7 ranges over all possible integer values of Xj. We repeat this until 
all subproblems are shown to be empty, or we find an integral point in one 
of them. 

B&B (and its version used to solve optimization problems) enhanced by 
cutting planes is a dependable algorithm implemented in most commercial 
software packages. However, instances in [HI [HJ [131 El 13 S] show that it is 
theoretically inefficient: it can take an exponential number of subproblems to 
prove the infeasibility of simple knapsack problems. While B&B is inefficient 
in the worst case, Cornuejols et al. in [TU] developed useful computational 
tools to give an early estimate on the size of the B&B tree in practice. 

Since IP feasibility is NP-complete, one can ask for polynomiality of a 
solution method only in fixed dimension. All algorithms that achieve such 
complexity rely on advanced techniques. The algorithms of Lenstra [22] 
and Kannan \Ta\ first round the polyhedron (i.e. apply a transformation to 
make it have a spherical appearance), then use basis reduction to reduce the 
problem to a provably small number of smaller dimensional subproblems. 
On the subproblems the algorithms are applied recursively, e.g. rounding 
is done again. Generalized basis reduction, proposed by Lovasz and Scarf 
in [23] avoids rounding, but needs to solve a sequence of linear programs to 
create the subproblems. 

There is a simpler way to use basis reduction in integer programming: 
preprocessing (pQ) to create an instance with short and near orthogonal 
columns in the constraint matrix, then simply feeding it to an IP solver. 
We describe two such methods that were proposed recently. We assume 
that A is an integral matrix with m rows, and n columns, and the Wi and 
£i are integral vectors. 

The rangespace reformulation of Q proposed by Krishnamoorthy and 
Pataki in [T7] is 



Findy €Q R nZ n , with Q/j = { y\ I < I ) Uy < 1 , (2) 




where U is a unimodular matrix computed to make the columns of the 
constraint matrix a reduced basis of the generated lattice. 

The nullspace reformulation of Aardal, Hurkens, and Lenstra proposed 
in [2], and further studied in [1] is applicable, when the rows of A are linearly 
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independent, and w\ = l\. It is 

Findy G Qn H Z n ~ m , with Q N = { y \ i 2 - x < By < w 2 - x } , (3) 

where xq E Z n satisfies Axq = l\ , and the columns of B are a reduced basis 
of the lattice {x £Z n \Ax = 0}. 

We analyze the use of Lenstra-Lenstra-Lovasz (LLL) |21j . and recip- 
rocal Korkhine-Zolotarev (RKZ) reduced bases [18] in the reformulations, 
and use Korkhine-Zolotarev (KZ) reduced bases [15] . |16j in our compu- 
tational study. The definitions of these reducedness concepts are given 
in Section [2j When Qr is computed using LLL reduction, we call it the 
LLL-rangespace reformulation of Q, and abusing notation we also call ([2]) 
the LLL-rangespace reformulation of (H|). Similarly we talk about LLL- 
nullspace, RKZ-rangespace, and RKZ-nullspace reformulations. 

Example 1. The polyhedron 

207 < 41xi + 38x 2 < 217 

( 4 ) 

< x h x 2 < 10 

is shown on the first picture of Figure [7J It is long and thin, and defines 
an infeasible, and relatively difficult integer feasibility problem for B&B, as 
branching on either x\ or x 2 yields 6 subproblems. Lenstra's and Kannan's 
algorithms would first transform this polyhedron to make it more spherical; 
generalized basis reduction would solve a sequence of linear programs to find 
the direction xi + x 2 along which the polyhedron is thin. 

The LLL-rangespace reformulation is 

207 < -3xi + 8x 2 < 217 

< -xi - 10x 2 < 10 (5) 
< xi + llx 2 < 10 

shown on the second picture of Figured now branching on y 2 proves integer 
infeasibility. (A similar example was given in 

The reformulation methods are easier to describe, than, say Lenstra's 
algorithm, and are also successful in practice in solving several classes of hard 
integer programs: see 0[TJ[T7]. -^ or instance, the original formulations of 
the marketshare problems of Cornuejols and Dawande in [9] are notoriously 
difficult for commercial solvers, while the nullspace reformulations are much 
easier to solve as shown by Aardal et al in [T]. 

However, they seem difficult to analyze in general. Aardal and Lenstra 
m [31 S] studied knapsack problems with a nonnegativity constraint, and 
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Figure 1: The polyhedron of Example [T] before and after the reformulation 

the constraint vector a having a given decomposition a = Xp + r, with p 
and r integral vectors, and A an integer, large compared to ||p|| and || r ||. 
They proved a lower bound on the norm of the last vector in the nullspace 
reformulation, and argued that branching on the corresponding variable will 
create a small number of B&B nodes. Krishnamoorthy and Pataki in [T7] 
pointed out a gap in this argument, and showed that branching on the 
constraint px in Q (which creates a small number of subproblems, as A is 
large), is equivalent to branching on the last variable in Qr and Qn- 

A result one may hope for is proving polynomiality of B&B on the refor- 
mulations of ([I]) when the dimension is fixed. While this seems difficult, we 
give a different, and perhaps even more surprising complexity analysis. It is 
in the spirit of Furst and Kannan's work in |12j on subset sum problems and 
builds on a generalization of their Lemma 1 to bound the fraction of integral 
matrices for which the shortest nonzero vectors of certain corresponding lat- 
tices are short. We also use an upper bound on the size of the B&B tree, 
which depends on the norms of the Gram-Schmidt vectors of the constraint 
matrix. We introduce necessary notation, and state our results, then give a 
comparison with [12J. 

When a statement is true for all, but at most a fraction of l/2 n of the 
elements of a set S, we say that it is true for almost all elements. The value 
of n will be clear from the context. Reverse B&B is B&B branching on the 
variables in reverse order starting with the one of highest index. We assume 
u>2 > ^2) and for simplicity of stating the results we also assume n > 5. For 
positive integers m, n and M we denote by G m ^ n (M) the set of matrices with 
m rows, and n columns, and the entries drawn from {1, . . . , M}. We denote 
by G' mn (M) the subset of G mjn (M) consisting of matrices with linearly 
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independent rows, and let 

XmAM) = \q-my (6) 

It is shown by Bourgain et. al in [6] that Xm,m(M) (and therefore also 
Xm,n(M) for m < n) are of the order 1 — o(l). In this paper we will use 
x(m,n, M) > 1/2 for simplicity. 

(A 

For matrices (and vectors) A and B, we write (A;B) for I I. For 

an m by n integral matrix A with independent rows we write gcd(A) for 
the greatest common divisor of the m by m sub determinants of A. If B&B 
generates at most one node at each level of the tree, we say that it solves 
an integer feasibility problem at the rootnode. 

If Q is a polyhedron, and z is an integral vector, then the width of Q 
along z is 

width(z, Q) = max { (z, x) \ x E Q} — min { (z, x) | x G Q} . (7) 
The main results of the paper follow. 

Theorem 1. There are positive constants d\ < 2, and di < 12 suc/t i/iai 
i/ie following hold. 

(1) If 

M > (dm \\(w 1 ;w 2 ) - (4;4)||) n/m+1 , (8) 

i/ien for almost all A £ G m , n (M) reverse B&B solves the RKZ-rangespace 
reformulation of (CP at the rootnode. 

(2) If 

M>(d2(n-m)\\w 2 -£ 2 \\) n/m , (9) 

then for almost all A £ G' mn (M) reverse B&B solves the RKZ-nullspace 
reformulation of |i]) at the rootnode. □ 

The proofs also show that when M obeys the above bounds, then Q has 
at most one element for almost all A £ G m ^ n (M). When n/m is fixed, and 
the problems are binary, the magnitude of M required is a polynomial in n. 

Theorem 2. The conclusions of Theorem^ hold for the LLL-reformulations, 
if the bounds on M are 

( 2 («+4)/2 \\( Wl] w 2 ) - (h-42)\\) n/m+ \ and ( 2 ("~ m + 4 )/ 2 || W2 -£ 2 ||) n / m , 

respectively. □ 
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Furst and Kannan, based on Lagarias' and Odlyzko's |19j and Frieze's 
work show that the subset sum problem is solvable in polynomial time 
using a simple iterative method for almost all weight vectors in { 1, ... , M} n , 
and all right hand sides, when M is sufficiently large, and a reduced basis 
of the orthogonal lattice of the weight vector is available. The lower bound 
on M is 2 cnlogn , when the basis is RKZ reduced, and 2 dn , when it is LLL 
reduced. Here c and d are positive constants. 

Theorems [T] and [2] generalize the solvability results from subset sum 
problems to bounded integer programs; also, we prove them via branch- 
and-bound, an algorithm considered inefficient from the theoretical point of 
view. 

Proposition Q] gives another indication why the reformulations are rel- 
atively easy. One can observe that &et{AA T ) can be quite large even for 
moderate values of M, if A E G m ^ n (M) is a random matrix with m < n, 
although we could not find any theoretical studies on the subject. For in- 
stance, for a random A 6 G4 3o(100) we found det(AA T ) to be of the order 
10 18 . ' 

While we cannot give a tight upper bound on the size of the B&B tree 
in terms of this determinant, we are able to bound the width of the refor- 
mulations along the last unit vector for any A (i.e. not just almost all). 

Proposition 1. If Qr and Qm are computed using RKZ reduction, then 

-,,w n \ s Vn \\(w 1 ;w 2 ) - nn . 
W1 dth(e n , Q R ) < det{AAT + I)1/(2n) ■ (10) 

Also, if A has independent rows, then 



n \ , gcd(A)y/n -m\\w 2 -£2 II n ,, 
w!dth(e n _ m , Q N ) < det{AAT)1/{2{n . m)) ■ (11) 



The same results hold for the LLL-reformulations, if ^Jn and i/n — m 
are replaced by 2( n_1 )/ 4 and 2( n_m_1 )/ 4 , respectively. □ 

Remark 1. As described in Section 5 of | '11], and in ]25^ for the nullspace 
reformulation, directions achieving the same widths exist in Q, and they 
can be quickly computed. For instance, if p is the last row ofU" 1 , then 
width(e n , Qr) = width(p, Q). 

A practitioner of integer programming may ask for the value of Theorems 
CD and EJ Proposition [2] and a computational study put these results into a 
more practical perspective. Proposition [2] shows that when m and n are not 
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too large, already fairly small values of M guarantee that the RKZ-nullspace 
reformulation (which has the smallest bound on M) of the majority of binary 
integer programs get solved at the rootnode. 

Proposition 2. Suppose that m and n are chosen according to Tabled and 
M is as shown in the third column. Then for at least 90% of A £ G' mn (M), 



n 


m 


M for 90 % 


M for 99 % 


30 


20 


33 


37 


50 


20 


1912 


2145 


50 


30 


96 


103 


60 


30 


420 


454 


70 


40 


197 


209 



Table 1: Values of M to make sure that the RKZ-nullspace reformulation 
of 90 or 99 % of the instances of type (fT2|) solve at the rootnode 

and all b right hand sides, reverse B&B solves the RKZ-nullspace reformu- 
lation of 

AX = b (12) 

x e {o,i} n 

at the rootnode. The same is true for 99% of A € G' mn (M), if M is as 
shown in the fourth column. □ 



Note that 2™~ m is the best upper bound one can give on the number 
of nodes when B&B is run on the original formulation fjl2[> : also, randomly 
generated IPs with n — m = 30 are nontrivial even for commercial solvers. 

According to Theorems Q] and [2l random integer programs with coeffi- 
cients drawn from { 1, . . . , M } should get easier, as M grows. Our com- 
putational study confirms this somewhat counterintuitive hypothesis on the 
family of marketshare problems of Cornuejols and Dawande in [9]. 

We generated twelve 5 by 40 matrices with entries drawn from {1, . . . , M} 
with M = 100, 1000, and 10000 (this is 36 matrices overall), set b = [Ae/2\ , 
where e is the vector of all ones, and constructed the instances of type (I12D . 
and 

b-e< Ax <b 
x£{0,l} n . 

The latter of these are a relaxed version, which correspond to trying to find 
an almost-equal market split. 
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Table [2] shows the average number of nodes that the commercial IP solver 
CPLEX 9.0 took to solve the rangespace reformulation of the inequality- and 
the nullspace reformulation of the equality constrained problems. 



M 


EQUALITY 


INEQUALITY 


100 


17531.92 


38884.92 


1000 


1254.42 


22899.67 


10000 


200.83 


1975.67 



Table 2: Average number of B&B nodes to solve the inequality- and 
equality-constrained marketshare problems 

Since RKZ reformulation is not implemented in any software that we 
know of, we used the Korkhine-Zolotarev (KZ) reduction routine from the 
NTL library [26J. For brevity we only report the number of B&B nodes, 
and not the actual computing times. 

All equality constrained instances turned out to be infeasible, except two, 
corresponding to M = 100. Among the inequality constrained problems 
there were fifteen feasible ones: all twelve with M = 100, and three with 
M = 1000. Since infeasible problems tend to be harder, this explains the 
more moderate decrease in difficulty as we go from M = 100 to M = 1000. 

Table [2] confirms the theoretical findings of the paper: the reformulations 
of random integer programs become easier as the size of the coefficients grow. 

In Section [2] we introduce further necessary notation, and give the proof 
of Theorems [T] and [2j 

2 Further notation, and proofs 

A lattice is a set of the form 

L = L(B) = {Bx\x e Z r }, (14) 

where B is a real matrix with r independent columns, called a basis of L, 
and r is called the rank of L. 

The euclidean norm of a shortest nonzero vector in L is denoted by 
Ai(L), and Hermite's constant is 

Cj = sup {Ai(L) 2 /(det L) 2/j \ L is a lattice of rank j ) . (15) 
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We define 

7 i = max{Ci,...,C i }. (16) 
A matrix A defines two lattices that we are interested in: 

L R (A) = L(A; /), L N {A) = {x e V l \Ax = 0} , (17) 

where we recall that (A; I) is the matrix obtained by stacking A on top of 
I. 

Given independent vectors b± , . . . , b r , the vectors b\,...,b* form the 
Gram-Schmidt orthogonalization of b\, . . . , b r , if b\ = b\, and b* is the 
projection of hi onto the orthogonal complement of the subspace spanned 
by bi, . . . , bi-i for i > 2. We have 

i-\ 

^ = b* ( 18 ) 

with 

Hi = (bi,b))/ || b* || 2 (l<j<i<r). (19) 
We call bi, . . . ,b r LLL-reduced if 

< ^(l<i<i<r), (20) 
> | || 2 (K*<r). (21) 
An LLL-reduced basis can be computed in polynomial time for varying n. 
Let 

i-l 

bi(k) =b* + Y, Vijbj (l<k<i<r), (22) 

j=k 

and for i = 1, . . . , r let Lj be the lattice generated by 

6i(i),6 i+ i («),..., 6 r (i). 

We call b\, . . . ,b r Korkhine-Zolotarev reduced (KZ-reduced for short) if bi(i) 
is the shortest nonzero vector in Lj for all i. Since L\ = L and &i(l) = b±, 
in a KZ-reduced basis the first vector is the shortest nonzero vector of L. 
Computing the shortest nonzero vector in a lattice is expected to be hard, 
though it is not known to be NP-hard. It can be done in polynomial time 
when the dimension is fixed, and so can be computing a KZ reduced basis. 

Given a lattice L its reciprocal lattice V is defined as 

L' = {ze lmL \ (z,x) € ZVx G L }. 



1/%' 

Hi,i-ibi-i + b * IP 
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For a basis b\,...,b r of the lattice L, there is a unique basis b'i, . . . , b' r of L' 
called the reciprocal basis of b\ , . . . , b r , with 



n & / \ _ 1 if i + j = r + l ^ 
otherwise. 



We call a basis bi,...,b r a reciprocal Korkhine Zolotarev (RKZ) basis of 
L, if its reciprocal basis is a KZ reduced basis of L' . Below we collect the 
important properties of RKZ and LLL reduced bases. 

Lemma 1. Suppose that bi,...,b r is a basis of the lattice L with Gram- 
Schmidt orthogonalization b\ , . . . , b* . Then 

(1) if b\, . . . ,b r is RKZ reduced, then 

||&?||>Ai(L)/Ci, (24) 

and 

\\b*\\> (detL) 1 ^'/^. (25) 

(2) if b\, . . . ,b r is LLL reduced, then 

II b* \\> Ai(L)/2^ -1 ^ 2 , (26) 

and 

\\b*\\> (detL) 1 / p /2( p - 1 )/ 4 . (27) 



Proof Statement (I24|) is proven in [18]. Let b±,...,b' r be the reciprocal 
basis. Since b[ is the shortest nonzero vector of L' , Minkowski's theorem 
implies 

||b'ill< v^(detL') 1/r . (28) 

Combining this with || b[ ||= 1/ || b* ||, and detL' = 1/detL prove (f25|) . 
Statement (]26p was proven in [21] . Multiplying the inequalities 

(29) 
□ 

(30) 



II b* || < 2( M )/ 2 1| b; || (i = i,... 

and using \\b*\\ ... || 6* ||= det L gives (p7|) . 
Lemma 2. Let P be a polyhedron 

P = {y eW \£ < By <w} 



and bi, . . . ,b* the Gram-Schmidt orthogonalization of the columns of B. 
When reverse B&B is applied to P, the number of nodes on the level of 
yi is at most 



( 


\\w-l\\ 


-)■■■( 


\\w-£\\ 




I WW \ 


WKW 



+ 1 



(31) 



10 



Proof First we show 



width(e r ,P) <\\w-£\\ / \\b*\\ . (32) 

Let x r) \ and x r ^. denote the maximum and the minimum of x r over P. 
Writing B for the matrix composed of the first r — 1 columns of B, and b r 
for the last column, it holds that there is xi,X2 £ W~ x such that Bx\+b r x r ^ 
and Bx2 + b r x r> 2 are in P. So 

1^ — ^11 > || (Bxi + b r x T ^i) — (Bx2 + b T x r ^2) || = || B(x± — X2) + b r {x r ^\ — x r ^) \\ 
> || b* || |x rj i — x T ^\ = || b* || width(e r , P) 

holds, and so does (1321), 



After branching on e r , . . . , ej+i, each subproblem is defined by a matrix 
formed of the first i columns of B, and bound vectors ti and Wi, which are 
translates of £ and w by the same vector. Hence the above proof implies 
that the width along ei in each of these subproblems is at most 

\\w-t\\/\\btl (33) 

and this completes the proof. □ 

Our Lemma [3] builds on Furst and Kannan's Lemma 1 in [12J, with 
inequality (|35|) also being a direct generalization. 

Lemma 3. For a positive integer k, let en be the fraction of A £ G mj „(M) 
with \i(Lr(A)) < k, andeN be the fraction of A G G' mn (M) with \i(Ln(A)) < 
k. Then 

(2k + l) n+m , . 

^ < . (34) 

and 

{2k + l) n 

CN - Mm XmA MY (35) 



Proof We first prove (I35p . For v, a fixed nonzero vector in Z n , consider 
the equation 

Av = 0. (36) 

There are at most M m ( ra_1 ) matrices in G' mn (M) that satisfy (|36|) : if the 
components of n — 1 columns of A are fixed, then the components of the 
column corresponding to a nonzero entry of v are determined from (|36|) . 
The number of vectors v in 7L n with || u ||< k is at most (2k + l) n , and the 
number of matrices in G' m n (M) is M mn Xm,n 

(M). Therefore 
„ (2k + l) n M m( - n ~V (2k + l) n 

£N ~ M mn X m,n( M ) ~ M m Xm ,n{ M )' 
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For ([3lj) . note that (v%) V2) € J J m + n \ s a nonzero vector in L R (A), iff V2 / 0, 
and 

Av 2 = v l . (37) 

An argument like the one in the proof of (|35[) shows that for fixed (v\; v 2 ) £ 
Z m+n with v 2 ^ 0, there are at most M m( - n ~^ matrices in G m , n (M) that 
satisfy (|37|) . The number of vectors in J J n + rn with norm at most k is at most 
(2k + l) n+m , so 

(2k + l)»+ m M m (" 4 ) (2k + l) n+m 



M mn M m 

□ 

Proof of Theorems [TJ and [2] For part (HJ) in Theorem [Q let 6J, . . . , 6* 

be the Gram-Schmidt orthogonalization of the columns of (A; I)U . Lemma 
[2] implies that the number of nodes generated by reverse B&B applied to 
Q R is at most one, if 

||6|||>||(^i ;U ; 2 )-(4;4)|| (38) 

for i = 1, . . . , n. Since the columns of (^4; I)U form an RKZ reduced basis 
of L R (A), dMD implies 

\\b*\\>Xi(L R (A))/Ci, (39) 

so (|38l) holds, when 

\ 1 (L R (A))>C i ||K;w2)-(*i;£ 2 )|| (40) 

does for i = 1, . . . , n, which is in turn implied by 

\ l (L R (A)) >ln ||K;to 2 )-(4;4)|| • (41) 

By Lemma l3l (pjTT) is true for all, but at most a fraction of e R of A £ G mjTl (M) 
if 

M (L27n||Wi;^2)-(^i;^)|| +lj)( m +^ 

l/m ' ^ ' 

and using the known estimate 7„ < 1 + re/4 (see for instance [24J) , setting 
e R = 1/2™, and doing some algebra yields the required result. 

The proof of part ([2]) of Theorem [T] is along the same lines: now 
&*,... ,bn- m is the Gram-Schmidt orthogonalization of the columns of B, 
which is an RKZ reduced basis of L^(A). Lemma [21 and the reducedness 
of B implies that the number of nodes generated by reverse B&B applied to 
Qn is at most one, if 

\ 1 (L N (A))>-/ n -. m \\w2-£ 2 \\, (43) 
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and by Lemma [3] (jl3j) is true for all, but at most a fraction of of A £ 
Gm,«(M) if 

M > (L^ II^-^II+U)^. (44) 

Then simple algebra, and using Xm,n{M) > 1/2 completes the proof. 

The proof of Theorem [2] is an almost verbatim copy, now using the 
estimate (|26l) to lower bound || b* ||. □ 



Proof of Proposition [Tj To see (fTUJ) , we start with 

width(e n ,P) <||w-*|| / ||6* || (45) 



from ()32p . combine it with the lower bound on || b* n || from (I25p . and the fact 
that 

det = det{AA T + I) 1 / 2 , (46) 

which follows from the definition of Lr(A). The proof of (|lip is analogous, 
but now we need to use 

det L R {A) = &et{AA T ) l l 2 / gcd(A), (47) 

whose proof can be found in [7J for instance. To prove the claims about the 
LLL-reformulations, we need to use ([271) in place of (|25l) . □ 



m r t%e^n°climenFionm^ ^n^ne^feWoliWp^o^s 

as an upper bound for N(n, k). The proof of Part ([2]) of Theorem factually 

implies that when 

M > M. 2 - m ?l™ 

£jV Xm,n{ M ) 

then for all, but at most a fraction of ejy of A G G min (M) reverse B&B 
solves the nullspace reformulation of (|12p at the rootnode. 



We use Xm,n{M) > 1/2, Blichfeldt's upper bound 

a<M^f. (49) 
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from [5] to bound Jn-m in (|48|) . dynamic programming to exactly find the 
values of N(n,k), and the values ejy = 0.1, and = 0.01 to obtain Table 

m 
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We note that in general N(n, k) is hard to compute, or find good upper 
bounds for; however for small values of n and k a simple dynamic program- 
ming algorithm finds the exact value quickly. □ 



Acknowledgement We thank Van Vu for pointing out reference [6]. 
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